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Abstract 

General formulas of the two-electron operator representing either atomic or effective interactions 
are given in a coupled tensorial form in relativistic approximation. The alternatives of using uncoupled, 
coupled and antisymmetric two-electron wave functions in constructing coupled tensorial form of the 
operator are studied. The second quantization technique is used. The considered operator acts in the 
space of states of open-subshell atoms. 

1 Introduction 

In the atomic structure calculations, investigations to optimize the effort of obtaining matrix elements 
of a two-electron operator are urgently required. This can be explained by the fact that theoretical 
methods recently used to produce high-precision atomic structure data generate large sets of matrix 
elements for a two-particle operator. This leads to large computation requirements in terms of both 
memory and speed. For example, large-scale configuration interaction (CI) calculations [TJ-[5] use 
a massive matrix for the atomic Hamiltonian. A large fraction of expansion terms of perturbation 
theory (PT) recently applied in atomic calculations are considered as matrix elements of some 

effective two-particle operators with complex tensorial structures. In all the above mentioned studies, 
a significant fraction of computations are devoted to the calculation of iV-electron angular parts of 
the matrix elements of a two-electron operator. Particularly complex calculations emerge when more 
than two open subshells with N > 2 are involved. A number of methods and techniques [12]- [17] were 
developed in order to obtain the general formulas for matrix elements of a two-particle operator for 
many electron case. The comprehensive description of this subject can be found everywhere [12] . |19| . 

In the present paper we distinguish the second quantization representation (SQR) [8], [13], [J3]. The 
efficiency of this technique manifests itself when the tensorial properties of creation and annihilation 
operators are taken into account [H]. Then the iV-electron angular part of a matrix element is 
described by a coupled tensorial product of creation and annihilation operators. In order to optimize 
(minimize) the calculation procedures it is important to choose the appropriate coupling schemes of 
angular momenta and the order of creation and annihilation operators in the tensorial product. In [5] 
and [16] the coupling schemes for the tensorial products of creation and annihilation operators were 
considered for nonrelativistic (LS'-coupling) and relativistic (jj-coupling) cases, respectively. The 
manner to determinate the expressions for matrix elements was presented. In [17] a coupled tensorial 
form of an effective two-particle operator used in a second-order MBPT was obtained in LS-coupling. 
Here, the different forms [TBj of coupling schemes to make tensorial product for particular cases were 
suggested. This enables one to reduce the complexity of the expressions for matrix elements. In [18] 
the investigations of [17] were extended by including into the presentation of coupled tensorial form of 
a two-particle operator coupled and antisymmetric two-electron wave functions given in LS-coupling. 
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In the present manuscript we continue the studies of [18] by considering a two-particle operator 
in jj-coupling (in the relativistic approximation [3]). We search for the general expressions of formal 
(effective) two-particle operator which describes atomic interactions as well as effective interactions in 
atoms. In Section [2] a coupled form of the two-electron operator is studied using uncoupled, coupled 
and antisymmetric two-electron wave functions. Sets of expressions for the two-electron operator are 
given in SQR (see Tables). An example of the application of obtained results for specific cases is 
considered in Section [3l 



2 Coupling schemes of ranks for a two-particle operator 

Let us consider a two-particle operator G given in the second quantization representation (SQR) in 
jj-coupling |14j : 

G = X ^2 aadpalal g(a,P,ti,is;-f,m y ). (1) 

In our considerations a, (3, u, fi indicate the subshells n a X a m a ,npXpmp, n u X u m u ,n fM X^m^ of iV-electron 

aX^ a A a n} ) X^ b Af ) ...nkX k k Af. (A ;,...) AM\\, the operator G acts on. Opera- 



wave function \^ N ) = 



tors aj and a\ denote electron creation and annihilation operators in the state niXiirii (Aj = liji, 
/j— parity of the state) with principal quantum number rij and magnetic quantum number (projec- 
tion) m-j. In the present paper the factor g(a, /3, fi, v\ 7, m 7 ) is associated with a matrix element 



n a X a m a nisX/3'm0 



5^(1,2) 



ni L Xi l mnn l/ X^'m u ) of a two-electron operator 



A' 



G 



i<j 



1 



2 9ij 



(2) 



of atomic interactions. We assume that can be expressed by the tensorial product [12 



9ij 
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(7) 



(3) 



71 72 



and that gij = gji. In ([3]) g(ri,rj) is the radial part of <?y. An irreducible tensorial operator g^\i) 
acts on the spin-angular variables of the i-th. electron in the space of one-electron relativistic wave 
functions (4-spinors) [19] 



\nXm) = \nljm) 



f(nlj\r) \ljm) 
-1) g(nl'j\r) \l'jm) 



(4) 



where I' = 2j — I, •& = I — j + 1/2. Functions f(nlj\r) and g(nl'j\r) are large and small compo- 
nents of \nXm). Functions \ljm), \l'jm) are 2-spinors. In (TjQ) the factor X = 1/2. However, when 
g(a, /3, fi, v\ 7, m 7 ) is associated with the antisymmetric matrix element |10] 



(5) 



(1,2) n fM X^m fl n v X v m u ) 
[1 - (/x o v)) (n a X a m a npXpm^ 5^ (1,2) n^X^m^n u X u m^ 



the factor X = 1/4. Here (/u <(->• v) indicates that [i must be interchanged with v. 

Notice, that when examining the atomic perturbation theory expansion terms or coupled clus- 
ter (CC) approach equation ones, they can be considered as the matrix elements of some effective 
operator e ^ g^} [6]. Then the factor g(a, /3, ju, u; 7, m 7 ) can be associated with the matrix element 



a/3 



grri~. 



[iv ) . The operator e ? f gml usually has more complicated tensorial structure and symme- 



try properties than g. 



(7) 



). Nevertheless, the expressions of G developed in this manuscript are also 



valid for e ^gml ■ In the later case the factor X is obtained individually. 
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Figure 1: The block-scheme of transforming two-electron operator into coupled tensorial form. 



Below we briefly describe the procedures that we have used to convert operator G into a coupled 
tensorial form, i.e., to obtain the expressions for G, where the quantities entering these expressions 
are independent of magnetic quantum numbers mj. It is convinient to explain such transformation by 
examining the schema (Figure [TJ which arises when applying a graphical method of angular momentum 
theory [21]. In schematic form we can write 



where summation indices a/3fj,u denote quantum numbers nX. The diagram A\ denotes either a 



matrix element ( a/3 



_(7) 
ym 1 



jiv ) or 



a/3 



ym 7 



fiv j. The irreducible tensorial product (the diagram 

A2) composed of creation and annihilation operators was produced by using Jucys theorems of graph- 
ical angular momentum theory |21j . To obtain the desired form of the irreducible product, we carried 
out several recouplings of angular momenta and made several changes of positions of creation and 
annihilation operators in (ff]). The arrangements of operators will be discussed later. Firstly, we shall 
discuss the recoupling of angular momenta. In this manuscript we investigate two approaches. In 
each approach the specific coupling schemes of the angular momenta j were applied. The block E of 
diagrams A3, A4 represents such schemes. In the hrst approach we have used the following coupling 
of momenta 



rn 



i< 



u 



3a 



Jv 



d 



while in the second approach, the block E is given by 



E K 



3a 



IP 



u 



J 11 
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Here [...] brackets define Clebsch-Gordan coefficients; u and d are the intermediate momenta arising 
in the recoupling procedure. The block T in A%, A3 describes the coupling schema of the irreducible 
tensorial product of creation and annihilation operators with the intermediate ranks Ji, J2. The 
diagram ^3 represents the recoupling coefficient of angular momenta ji transforming the schema E 
into T. In the case of ([3]), due to the chosen specific coupling £5, in the first approach the diagram A4 
corresponds to product of the submatrix element [A Q A/3 \\g^ || A^A^] and 6(u, jx)S(d, 72). However, in 
the second approach A^ is associated with submatrix element 



ff W(l,2; 



A/jAjyd 



J2 [ja,jp,J,d} 1/2 
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A a A/3 



,(7172) 



A/jAz, 



Rafiiivi)-, 2). 

(9) 
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Here i? Q/ g /i;/ (l, 2) is a radial integral of a radial function g{r 1^2) in the basis of \nlj) functions. In 
([9]) the coupled two-electron wave functions 



\rii\inj\jum u ) = 



Aj Xj u 



\mXiminjXjmXj) (10) 



are used to determine the matrix element of gml • Note that in the first approach ([7]) uncoupled wave 
functions \n a X a ra a n^X^rap) = \n a X a m a ) \npXpmp) were employed. Bellow the indices b and z denote 
the quantities which have been obtained in the first ([7]) and the second ([8]) approaches, respectively. 

The coupling scheme of the block T is developed to consider the order of creation and annihilation 
operators (A2). Let us study this problem in a more detail. We collect the terms of operator G taking 
into account on how many subshells of equivalent electrons creation and annihilation operators act on. 
Then we can write 

G = + ^Cjj + Gijk + Gijki- (11) 

i i<j i<j<k i<j<k<l 

Operators Gi, Gij, Gijk, G^ki act in the space of the states of one, two, three and four subshells, 
respectively. Indices k, I numerate the subshells in \^ N } the operator G acts on. In our study, the 
sums in (jllh run in a way that i < j < k < I. The placing (arrangement) of creation and annihilation 
operators in (fTTj) follows the suggestions of |17j : first of all, operators ami an d a>m{ which act on the 
same subshell are collected side by side; secondly, operators a,m\ and cim\ acting on the first (second, 
third) subshell of many-electron wave function are situated to the left of the ones acting on the second 
(third, fourth) subshell. Each operator in (llip is given in the coupled form 

G L .,= Y, *G& J *^*g{a,Q e ,J u J 2 ,i). (12) 

Jl J2 Be 

x G^ J2 ^ 7 m 7 denotes the irreducible tensorial product of the operators Qm\ and am\ with the in- 
termediate ranks J%, J 2 and with the resulting rank 7. T (see diagram .A3) defines the coupling 
scheme of the irreducible tensorial product. The factor x g(s, g e , Ji, J2, 7) (the diagrams A3 and 
A4) includes the submatrix elements of g^ and the recoupling coefficients arising while making 
the tensorial product. In (|12p superscript x prescribes in which approach the quantities are ob- 
tained. Argument x = b, z indicates the first and second approaches, respectively; s and g e char- 
acterize the set {n a X a n^Xpn^X^n v X v } of quantum numbers in (fT2]k s indicates the number of sub- 
shells the operator G acts on. We collect together the terms of operator G with definite s into 
the groups. Operators G^ lJ2 ^ 7 m 7 of a particular group connect exclusively the configuration states 
(. . . JVi, . . . N2, ■ ■ ■ N$, ...N4 . . .\ and |. . . JV(, . . . N' 2l . . . ...JV4 . . .) with the specific electron occupation 

numbers iV, and iV-, i.e., {5 1 = N{-N 1 , 5 2 = N^-N 2 , 5 3 = N^-N 3 , 5 A = N4—N4}. Furthermore, the 
terms of each group g are collected into the subgroups (e numerates different subgroups). Each sub- 
group g e is characterized by the following sets of the quantum numbers: stl = {noAQn^A^n^A^n^A^}, 
st2 = {npXpriaXan^X^nyXy} , stS = {naXaiipXpn^X^nyXy} , sti = {n^A^n^Aan^A^n^A^}. The sub- 
groups differ the sets with distinct collections of n a X a , npXp,n u X u ,n^X^ quantum numbers. Note 
that the terms with stl (st2) and st3 (stA) describe the direct and exchange interactions, correspond- 
ingly. Creation and annihilation operators with fixed g e compose the irreducible tensorial product 
G^sg e ^ 7 m 7 (the diagram A2). In general, the factor x g(s, g e , J\, J2, 7) has four terms associated 
with stl, st2, st3, st4. However, due to the symmetry properties of atomic interactions, the term in 
x g(s, g e , Ji, J2,7) corresponding to the set stl (st3) is equal to the term described by st2(st4). 

For convenience, the expressions of x G^sq/ 2 m 7 and x g(s, g e , Ji, J2,l) for s = 2, 3, 4 are collected 
in Tables [DCil The expressions for the operator G when it acts on one subshell (n a X a = npXp = 
n^X^ = n v X v , s = 1) can be found everywhere (see, for example, in [2], [16], [17]), therefore they are 
not presented in this paper. 

Let us concentrate now on the study of the tensorial part x G^q/ m of G for s > 1. De- 
creasing a number of expressions which should be written for G^ m , we explore the convention 



4 



Table 1: The quantities for generation of the expressions for the operator G in two-subshell case. 
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(0,0) 
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1 


T 21 (-l, 1,-1,1) 


x((-l)^6i 2 i 2 ,D + 0i) 
+ (b 1221 D b 2 + fa) 


— ^1212-021 


(-2,2) 


2 


1 


T 21 (-l, -1,1,1) 




■^1122 -D 22 


(2,-2) 


5 


1 


T 21 (l, 1,-1,-1) 


6 5 , 22i 


2 5'22i 


(-1,1) 


3 


1 


*T 22 (-1, -1,1,1) 


(_l)J2+J2+7i+72('5 

x£j l72J P + 0l) 


^1112^Jim-D23 






o 
Z 


^T 1 ('I 1 1 1 ^ 
J 23 1 — 1 , — -L, J- J 


(-1) J1+J2+7 (&1222 

x5j 17a P|j m _,- a + 4>i) 


-Zl222<5jid 


(1,-1) 


4 


1 


K T 22 (-1,1,1,-1)| 6 ^ 


(-1)*-* b g 23l 


(_l),W2^ 23i 


2 


^^(l, — 1, — 1, 1) 1 6-f-*-* 


(-1)*-* 0<7 232 


(_l)jw a ^ 23a 



T a p(bi, b 2 , 63, 64) to describe the tensorial products of creation and annihilation operators. In T a p the 
i-th operator is equal to or if the argument bi = —1 or 1. For instance, T 2 i(&i, 6 2 , 63, 64) 
describes three operators G^g 



'"" ' J ~ " 1 when g = 1,2 and 5 (see Table [T]). It is important to note that 



in our study in the case of three and four subshells the operators x G^ g 1 ^ 2 ^ m are identical for x = b 

fe^<(JiJ2)(7) anc j zq(JiJ2){i) 



2g e 



differ. More 



and x = z. Only for two-subshell case, the operators G 2Se 

exactly, this takes place only if {61,62} = {—1, 1}, {1, —1}, (i.e., when g = 3,4) and then e = 1,2 for 
ixed g. Thus, for 
x g(s,g e , Jl, J 2 ,7). 

Operators G 2 ^ 
products with the following coupling schemes: 

(Ai.Ai) xT^I (a 2 ,A 2 



fixed £. Thus, for the remaining cases index e is redundant and it is dropped out in x G^f 2 ^ 1 m 7 and 
Operators G^ g x f m ^ (Table [[]) which act on two subshells (s = 2) are described by tensorial 



T 21 (-l, 1,-1,1) := 



6 T 22 (-1, -1,1,1) : 
6 T 23 (-1, -1,1,1) : 
Z T 22 (-1, -1,1,1) 



(7) 



(J2) 

(a 2 ,A 2 ) x a 

(■&) 



«WxrW fAi, Ai 



x 



x a 



(As) 



*T 23 (-1, -1,1,1) : 
Here we have used the definition 

T« (A, A,): 



T( Jl ) (Ai,Ai) xS^ 

a (A 2 ) xT W(A 2 ,A 2 



x a 



x 



/ J m 


7 


i(A 2 )" 


(7) 




rrij 


wo - 


(7) 




m 7 


(A2) 


(7) 






MO" 


(7) 




m 7 



x a^) 



M) 



(13) 
(14) 
(15) 
(16) 
(17) 

(18) 



For three subshells (s = 3), operators G 



tensorial products: 



(JiJb)(7) 
3g m-y 



(Table [2J are given by the following types of 



T 31 (-l, 1,1,-1) := 



(\ 1 ,X 1 ) 



(Ja) 



x a 



(A3) 



(7) 



(19) 
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Table 2: The quantities for generation of the expressions for the operator G in three-subshell case. 
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(0 -1 1) 
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Toif— 1 1 —1 11 

- L 31\ - L ,- L , x , -v 


-(6 1213 £iV + 1 ) 

+ (&1231<W TiV + 2 ) 


— ■^1213 
X -D 3 i 


(0,1,-1) 


1 


T 3 i(l, -1,1,-1) 


("l) Jl 6 <?32 


("l) Jl ^32 


(-2,1,1) 


4 


T 31 (-l, -1,1,1) 


(b 1123J D 3 b + 00 


■^1123-D 33 


(2,-1,-1) 


3 


T 31 (l, 1,-1,-1) 


^34 














(-1,1,0) 


5 


T 32 (-l, 1,-1,1) 


x((-l)^6 1323 PiV + </»!) 

+ (^1332^2 + 02) 


■Z'l323-D 3 2 


(1,-1,0) 


6 


T 32 (l, -1,-1,1) 


("1) J2 6 fe 


("1) J2 ^35 


(-1,-1,2) 


7 


T 32 (-l, -1,1,1) 


(b 1233 Piv + 00 


-^1233^22 


(1,1,-2) 


8 


r 32 (i, 1,-1,-1) 


6 <?37 


^37 


(\ n — il 


Q 


Too (4 —1 1 —11 


x 6 93il • 

^ 6 I.71-B-.72 


(_l)ii+^i-^2 
x z 93i\ h ^ j2 


l u , J-j 


i n 


Toof— 1 1 —1 11 


X ^032 1 • 

I.71-B-.72 




f 1 9 11 


1 1 
1 1 


J-33\ — -L, J-, -L, — J- J 


I.71-H-.72 


X ^Sslj^^ 


(1,-2,1) 


12 


T 33 (l, -1,-1,1) 


(—lyx+Ji-Ji 

x 6 g 34 | . 

yJ4 I.71-H-.72 


(_l)ii+^i-^2 
x *fl34| il++ia 



T 32 (-l, 1,-1,1) := M Jl ) (Ai,A 2 ) xT< J2 > (a 3 ,A 3 



(t) 



and 



T 33 (l, -1,1,-1) := 



^'xTW fA 2 ,A 2 



x a 



(As) 



(7) 



Finally, the four-subshell case (s = 4) is presented by the tensorial product 



T 41 (-l, -1,1,1) := 



(Ai,A 2 ) xa^ 



M0 



x a 



(7) 



(20) 
(21) 



(22) 



and operators C^ 1,72 ^ 7 ^ in Table Note again, that in (fT3"]) - ([2"2"j) operators T Q( g(6i, 6 2 , 6 3 , 64) de- 
scribe all types (in the sense of coupling scheme) of the tensorial products used in the present pa- 
per. The expressions for operators X G^ J2 ^ 7 m are easily obtained from ^3|) - ([22]) and Tables [Lj3j 
When a consecutive coupling of the resulting moments of subshells in many-electron wave function 

a X^ a A a nb\^ b Ab...nk\u k (A a fc...) AMa) is used, the formulas for matrix elements of x G^ s Jl " 1 1 



n 



can be immediately found from general expressions given in [17] . 

Consider now in detail the factor x g(s, g e , J\, J 2 , 7) (fl2j) . When x = b (the first approach), we 
obtain 

b g(s,Q e ,J 1 ,J 2 ,j) = ^^({ b D SQe ('y, X a , Xp, X v , A M 

, 7i, 72, J\i J2) 

7172 

x 6(7, A Q , Xfs, X„, A At ,7i,7 2 ) + </>i) + ( fc -E S£ , e (7,A a ,A ( g,A J ,,A / ,,7i,7 2 , Ji, J 2 ) 

x 6(7, A Q , A/?, A^, A^, 71,72) + </> 2 )). (23) 
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Table 3: The quantities for generation of the expressions for the operator G in four-subshell case. 
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+ 02) 
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(-1,1,1,-1) 
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T 41 (-l, 1,1,-1) 




^43 



The expressions for the special cases of ft <?(s, g e ,Ji, ^2,7) are given in the fifth column of Table [Hand 
the fourth column of Tables [21 [3j In the case of ([3]) 



Kpvn = Kt, ^aA/3,K, A M ,7i,7 2 ) := 





1/2 










A a A/3 


0(7172) 


A^A^ 


.71,72. 







i2a^ (1,2). (24) 



The factors 6 L> S g e and 6 S see arrive due to the recoupling procedures described previously. The 
explicit expressions of these factors are obtained by using the relations in Table H] dealing with the 
following ones: 



z 



X ~1y 

Ixy y 



and 



1 /2 

L(x 1 ,x 2 ,yi,y2,7x,7y,lxy,zi,z 2 ) := [72,72/^1^2] 7 



xi x 2 7x 



yi 

21 



V2 ly 
Z 2 Ixy 



(25) 



(26) 



To give more compact expression for b g(s, g e , Ji, J 2 , 7), we have introduced the notations 0i and 
2 for the second terms in the brackets in (I23D . The formula for the term 4>i (</> 2 ) linked to st2 
(si4) is obtained from the expression of the first term in the brackets associated with stl (st3) by 
replacing b a p vli (6^0 with bp aiiV (bpau^), interchanging 71 and 7 2 (71 o 72), and_multiplying the 
obtained formula by the factor (—1)71+72-7 = w _ Furthermore, in Tables the notation b g(s, g' e , J\, J 2 , 7) 
implies that the expression for b g(s, g e , J\, J 2 , 7) is found from one of b g(s, g' e , Ji, J 2 , 7) when b a $ Vil 
is replaced with b^ v p a and the obtained formula is multiplied by (— l)-? a +J/3+>+-?M-7i-72_ Note that 

2) is equal to the term associated with 

1, 

1/2 



for the atomic interactions when gu = gji, the term <fi] 

stl (st3). When some effective interaction e ^g[l m „ is studied, the expression for G can be found 



by replacing 6(7, A a , A/?, Aj,^^, 71, 7 2 ) with the factor ^£ [A^A^ \\ eSf g^ \\ \^K] ■ Finally, for 

obtaining the expression of b g(s, g e , J\ } J 2 , 7) in the case of the antisymmetric matrix element ([5]), the 
factor 6(7, X a , A/3, X u , A M , 71, 72) of the term with stl in ([23]) must be replaced by 



-(6(7, A a , A/3, A„, A M ,7i,7 2 ) 



Pl,P2 
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Table 4: The relations for determination of the recoupling coefficients. 
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Xi^O'l, J 2 ,72, il,j2, Jl) 


PN 


(-l) 1 - j2 L(ji,j3,j 2 ,i 3 ,7i,7 2 ,7, Ji, J 2 ) 


Four-subshell 






(_l) Jl - J 4-7l+7K( Jl , j4 , 72 , 7l , 7 , J 2 ) 
Xi^O'l, J 2 ,7l,j2,j3, A) 


U V2 


(-I) j3+j4 - 71 ^(j3,j4,7l,72,7,^2) 


D z 


5( k( J3) i4) d) u> 7) j 2 ) 


D z 

U V2 


(_l)-+«+7+I^ 1>J - 4>djUj7j J 2 ) 
Xi^O'i, J 2 iU,j 2 J 3 , Ji) 


DV 


(_l)-7+^4+ JlK(i72)i4)72)7l5 75 J 2 ) 

xX(j 2 , J 2 ,7i,ii,j3,^i) 
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3u 3a 71 

x{ jp. 3n 72 > b(j,X a ,Xp, X ll ,X u ,p l ,p 2 )). (27) 

Pi P2 7 J 

In addition, the terms corresponding to st2, st3 and stA must be dropped in ([23 p. 
In the second approach, the general formula for x <7 (x = z) is presented as follows: 

z g(s,g e ,J 1 ,J 2 ,'y) := ^ z E SQe (J 1 , J 2 , 7, A a , A^, \„, A^, u, d)Z(7, A a , A^, A„, A M , u, d). (28) 



The expressions of the factors z E sge (Ji, J2,7, X a ,Xp,X v , A^,u, d) which represent the recoupling co- 
efficients are given in the last column of Tables QJ3J The second factor on the right side of (|28p is 
expressed as 

Zafivii = X a ,Xp,Xu, Xfj,,u, d) := z a/3ufl + r](a, /3, u)z/3 aufl + r](^, v, d)z a/3liu + r](a, ft, u)t/(/x, v, d)zp a ^ v , 

(29) 

where rj(x,y,f) := (1 - 5(x,y))(-l) jx+jv ~~ f+1 . In the case of ©, 



^M^^Aa.A^A^A^.ti.d) := ^(-l)^+>- rf [n/ 7 ] 1/2 X a Xpu ^(1,2) 



XnXyd 



(30) 



Note that when = ^ the fourth (third) term in (j29p gives the same contribution as the 
first (second) one. For antisymmetric coupled two-electron wave functions, the submatrix element of 
operator </^ 7 is expressed as 



X^X v d 



[l_(_l)^+>-7 (At0z /)] l Xa Xpu 5 (7)(i ;2 ) 



X^X v d 



and 



1 



^(7, A a , A/?, A„, A^, it, d) := - (-1) J 



Furthemore, when studying the effective interaction 

z( 7 ,A a ,A^,A„A^,n,d) := [u/7 ]i/2 r VA/3U 



It 


1/2 










XaXpU 


5 (7) (1,2) 


X^Xi/d 


.7. 









XaXyd 



(31) 
(32) 

(33) 



In Tables the notation z g(ce, g' , Ji, J2,l) means that the expression for z g(a, g, J\, J2,l) is found from 
the expression for z g(a, g', J\, J2,l)\u-<^d, where z a p ull is replaced with z^pa and the obtained formula 
is multiplied by the phase factor (— l)-? a +.W+>+.?M-7 - The explicite expressions of (|24p . ()33|) submatrix 
elements of Coulomb, magnetic and retardation interactions can be found, for instance, in [12j . |19j . 



3 Closed subshell cases 

Particularly simple expressions for the two-electron operator ([1]) can be obtained when it acts on at 
least one closed subshell (say n a Xa a+1 ) of j 1 I /Ar ). Then, in the first approach from ([I]) we obtain 

h G afill = N a T^(Xp,X^)[j a ]- l Y[-[j a ] 1/2 (5 ll0 5 l2l b a p^ a + 5-y 1 'y5-- > -20^/3aa/i) 



7i 72 



+ ( _ 1) i, + i, + 7 [7l)72] l/ 2x (\ J J_ J \ bftaiia + u \ I J J ) b a ; 



3ii 3/3 3c 



3p 3n 3c 



(34) 



In the second approach 

z G a p„ = 2 Y, ^^(A^A^bal-^^lf+^^d] 1 ^/ u d ])z a p m . (35) 

x „ n I 3ft 3/3 3a ) 
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Here the operator TV = —[A] 1 / 2 [a^ x a^]^ has a submatrix element equal to N. In the case 
when G acts on two closed subshells (say n a \ 2 j a+1 , np\ 2 p l3+1 ), we obtain 

b G a p= Y, ^o^^LW/j] -172 ^ {<5 7 i0^ (36) 

npXp 7172 

z G a p = -2 Y, ^oNaNpljaJp^YSudM^ZaPap- (37) 

Finally, when G acts on one closed subshell (say n a X^ a+1 ), we obtain 

b G a = Y Y^ a ^ a ~ x ) ^70^71 0^720 [j a] baaaa- 

(38) 

n a X a 7172 

*G a = Y fi<*(N a - l)5 lQ [j a ]- 2 Y^ud{-l) u [u] l,2 z aaaa . (39) 

4 Conclusions 

A two-electron operator G which describes the relativistic interactions in atoms was considered in a 
coupled tensorial form in 7 7-coupling. The second-quantization representation was used. A complete 
set of the expressions when G acts on two, three and four subshells (the largest number of subshells 
the operator can act at the same time) of many-electron wave function \^ N ) are presented in a 
compact form. It allows easy generation of the formula for G when the particular case is considered. 
Each expression is given in such a structure that the calculation of the matrix elements of G can 
be performed into two separate tasks: the calculation of TV-electron spin-angular part presented by 
a submatrix element of the irreducible tensorial product G^g^^m-, composed from creation and 
annihilation operators (the most computer time-consuming part of calculations) and the determination 
of the factors x g(s, g e , J±, J2, 7). The factors x g(s, g e , J\, J2, 7) do not depend on the resulting angular 
momentum of subshells A a of | , I /Ar ), thus they can be determined before the calculations of submatrix 

elements for G^sg'^^m-, and reduce the computation time of matrix elements. 

In the present paper we apply the coupling schemes for G^ lJ2 ^ 7 m 7 which are very useful in the case 
of the consecutive order coupling of the resulting momenta of subshells in J^^. Then many-electron 

submatrix element of G ( / g lJ2 ^' 1 m 1 takes a very simple expression, i.e., it is expressed as sums which 
run over only intermediate ranks (J±, J2) of the products of 6j and/or 9j-symbols and submatrix 
elements of operators acting in the space of states formed from a subshell of equivalent electrons. 
In the present paper two forms of coupling schemes of angular momenta of two-electron operator G 
were studied. It enables us to use uncoupled (the first approach), coupled (the second approach) and 
antisymmetric two-electron wave functions in constructing coupled tensorial form of the operator. The 
possibility to apply different types of two-electron wave functions allows one to choose more optimal 
ways of calculations. Note that the second approach is preferable for the problems where several 
operators with different tensorial structures are considered, for instance, in the formation of energy 
matrix for the atomic Hamiltonian. In this case, Coulomb and Breit interactions can be presented 
by a single operator G with two-electron submatrix element \\ a \pu ||<7^|| A^A^d] when 7 = 0, where 
g(°) = g$ ou i omb + g^Breif Then many-electron angular part can be determined as the submatrix element 
of the unique operator G^sg^^m-, ■ The first approach is more preferable when one seeks to calculate 
the matrix elements of particular operator efficiently. In this approach, the internal tensorial structure 
of the operator g^) is directly involved (through the diagram A4, the coupling scheme Eb) into the 

expressions of many-electron matrix element of G^ lJ2 ^ 7 m 7 and recoupling coefficient (the diagram 
A A ). 

It is important to note that the expressions of both approaches are also applicable to the study of 
the operators representing some effective interactions in atoms arising, for instance, in Atomic MBPT 
or Coupled Cluster (CC) method. 
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The method to obtain the formulas for the operator G developed in [3J, [16], can be explained as the 
combination of our first and second approaches. However, the methodology proposed in our paper on 
coupling schemes to construct the irreducible tensorial products of creation and annihilation operators 
allows us to find expressions for many-electron matrix element of the operator G more simply than in 
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